Clustering and Triaxial Deformations of ^°Ca 
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We have studied the positive-parity states of '"^Ca using antisymmetrized molecular dynamics 
(AMD) and the generator coordinate method (GCM). Imposing two diflterent kinds of constraints 
on the variational calculation, we have found various kinds of ^''Ca structures such as a deformed- 
shell structure, as well as a- At and C- Si cluster structures. After the GCM calculation, 
we obtained a normal-deformed band and a superdeformed band together with their side bands 
associated with triaxial deformation. The calculated B{E2) values agreed well with empirical data. 
It was also found that the normal-deformed and superdeformed bands have a non-negligible Q:-"^^Ar 
cluster component and ^^C-^**Si cluster component, respectively. This leads to the presence of an 
a-^^Ar higher-nodal band occurring above the normal-deformed band. 

PACS numbers: 21.60.-n, 23.20.-g 



I. INTRODUCTION 

Nuclear dynamics possess various aspects depending 
on mass regions, excitation energies, and so on. In light- 
weight nuclei, it is known that clustering plays a signif- 
icant role in the features of ground and excited states 
[l|, On the other hand, in heavier nuclei, the cluster- 
ing effects are not clear, though many theoretical 0, 0| 
and experimental Q studies have been conducted. In 
the /p-shell region, the focus should be on proton-rich 
N Z nuclei, because such nuclei can have a clustered 
structure comprising stable nucleus. Moreover, proton- 
rich nuclei have a large radius and stronger Coulomb re- 
pulsion, and may also derive a cluster structure. In the 
scope of this research, the features of cluster structures 
in '*°Ca are a key issue, because they are the heaviest 
N = Z stable nuclei, and many experimental data ex- 
ist for this nuclei. In this paper, we have studied the 
structure of ^"Ca as a starting point to understand the 
structures of medium- and heavy-weight N ^ Z nuclei. 
^°Ca has a typical double closed-shell structure nucleus 
and has a spherical ground state. However, it is known 
that many kinds of deformed band appear in low energy 
regions. The first K"" = 0+ band built on the J" = 
state (3.35 MeV) is considered to be a normal-deformed 
(ND) state and the dominant configuration is 4p-4/i Q. 
The K'^ = 2+ band built on J'^ = 2+ (5.25 MeV) exists 
just above the K"^ = 0+ band. It has been suggested that 
the ND band deforms triaxially and has the K"^ = 2+ side 
band due to triaxiality 0, [1] . 

The a-'^^Ar cluster structure has been studied for a 
long time, because ^"Ca is an analogue of ^^O, which has 
an a-^^C cluster structure in the first K'^ — 0+ band, as 
a double closed-shell nuclei. The local potential model 
|JllJn| and a-^^Ar orthogonal condition model (OCM) 
12 ■ Il3l| have been performed theoretically. Ohkubo et 
al. suggested that the first K"^ = 0+ band (ND) has an 
a-^^Ar structure, and predicted that its parity-doublet 



= band and a-'^^Ar higher nodal band exist in 
highly excited states Sakuda et al. obtained the 

K"" = 2+ band as well as the K"" = 0+ and 0~ states 
using the a-'^^Ar OCM^ and succeeded in reproducing 



E2 transition strengths |13|. Experimentally, the a-'^^Ar 
structure is studied through the ^^Ar(^Li,d)*°Ca reac- 
tion [11 [H- The states in these K"" = 0+ and 0" bands 
are populated by the a-transfer reactions and have large 
a spectroscopic factors fl^. In the experiments we de- 
scribe, the a-'^^Ar higher-nodal states were also observed 

m- 

It has been suggested that the states in the K'^ = O"*" 
rotational band built on the J"^ = state (5.21 MeV) 
have a 8p-8h configuration 0], and have been observed 
during experimental work searc hing for the 8p-Sh states 
with ■^^S(^^C, a)^°Ca reactions [ll]. Due to the strong 
population in the multi-nucleon transfer data and the 
strong E2 transitions [TtI, the 0+ (5.21 MeV), 2+ (5.63 
MeV) and 4+ (6.54 MeV) bands have been thought to 
belong to the superdeformed (SD) band with the domi- 
nant 8p-8h configuration. Recently, by using GAMMA- 
SPHERE array detectors, the level structure of the de- 
formed bands in ^"Ca has been explored and many ex- 
cited states up to high spin have been discovered. This 
band was thus confirmed as the SD band p^ . 

Motivated by these the experimental observations, 
many theoretical microscopic studies on deformed states 
of ^"^Ca have been performed recently with the methods 
of Skyrme-Hartree-Fock (SHE) [l|, SHF-BCS + GCM 
[lOl, spherical-basis AMD i2l| and the shell model [H]. 
Inakura et al. performed cranked SHE calculations with- 
out assuming axial symmetry, though energy levels were 
not calculated JJJ. Bender et al. performed SHE-BCS -I- 
GCM calculations [l^l • Although they calculated energy 
levels and quadrupole transition strengths in the ND and 
SD bands, they could not study triaxiality nor side bands 
because they assumed axial symmetry. In these stud- 
ies, the relationship between deformed states and clus- 
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ter structure was not discussed. It has been suggested 
that the SD state forms a ^^C-^^Si-hke cluster structure 
in spherical-basis AMD [2l|. Within the spherical-basis 
AMD, triaxiality does not appear in ND nor SD states. 

The purpose of the present study is to understand the 
clustering and triaxial deformations in the low energy 
states of '^'^Ca in a unified manner. We use the frame- 
work of AMD + GCM. The basis functions of GCM are 
obtained by energy variation after parity projection with 
constraints. We adopted two kinds of constraints. One is 
a constraint on the quadrupole deformation parameter P 
(/3-constraint) and the other is on the distance d between 
clusters' centers of mass (d-constraint). It has already 
been proven that the d-constraint is useful for obtaining 
various kinds of clustering wave function, which arc not 
computed within a simple /3-constraint 23]. For exam- 
ple, the ^Be(2a)-^^C cluster structure in ^°Ne is calcu- 
lated with a d-constraint but not with the /3-constraint. 
Also in the case of ^°Ca, many kinds of cluster struc- 
ture, for example a-'^^Ar, ^Be(2a)-'^^S and ^^C-^*Si, can 
be calculated with the d-constraint, although no cluster 
structure is obtained in '"'Ca with the /3-constraint. We 
superposed mean-field-type and cluster-type wave func- 
tions calculated with /3- and d-constraints respectively, 
and calculated energies and E2 transition strength. We 
analyzed the superposed wave functions in order to in- 
vestigate clustering and triaxial deformations. 

This paper is organized as follows. In the next section 
fini|. we explain the framework of this study. The cal- 
culated results and discussions are presented in ijllll and 
lastly, we present a summary in ijlVI 



part \xi) is parametrized by a; and /3i and the isospin 
part \Ti) is fixed as \p) (proton) or |n) (neutron). The 
quantities (Z^, a^, /3i, ^'j;, ^'J^, z/^) are variational parame- 
ters and are optimized by energy variation as explained 
in the next subsection. 

The trial wave function in the energy variation with 
constraints is a parity-projected wave function. 



^int), 



(2) 



where tt is parity and Pr is the parity operator. In this 
study, we will discuss positive parity states. 
The Hamiltonian is. 



H = K + V^ + Vc-Kc 



(3) 



where K and Kq are the kinetic energy and the energy of 
the center of mass motion respectively, and Vn is the ef- 
fective nucleon-nucleon interaction. We have used Gogny 
DIS force (DIS) and Skyrme SLy7 force (SLy7) in the 
present work. The Coulomb force Vc is approximated by 
a sum of seven Gaussians. 



B. Energy Variation, Angular Momentum 
Projection and the Generator Coordinate Method 

We performed energy variation and optimized the vari- 
ational parameters included in the trial wave function 
(Eqs. H])) to find the state that minimizes the energy of 
the system E'^, 



II. FRAMEWORK 
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A. Wave Function and Hamiltonian 

We used the theoretical framework of AMD + GCM. 
In the present study, the AMD wave function is a Slater 
determinant of triaxially deformed Gaussian wave pack- 
ets (deformed-basis AMD), 
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Here, the complex parameters Z^, which represent the 
centroids of the Gaussian in phase space, take indepen- 
dent values for each single particle wave function. The 
width parameters Vx, and v^, are real parameters and 
take independent values for each of the x-, y- and z- 
directions, but are common for all nucleons. The spin 



Here, we add the constraint potential Vcnst to the expec- 
tation value of Hamiltonian H in order to obtain the min- 
imum energy state under the optional constraint condi- 
tion. In this study, we employed two types of constraint, 
which are on the quadrupole deformation parameter /3 
(/3-constraint) and the distance between clusters' cen- 
ters of mass, d (d-constraint) by employing the potential 

^nst 1 



"cnst 
^cnst 



rfc„,-c„)' 



for /3-constraint, 
for d-constraint. 



(5) 

Here /3 is the matter quadrupole deformation parameter, 
which is defined in Ref. [2J, and dc„-c„ is the distance 
between the clusters' centers of mass Cm and C„, 



IRn 



(6) 



(7) 



where A„ is the mass number of cluster C„ and the ex- 
pression i G Cn means that the zth nucleon is contained in 
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cluster C„. It should be noted that the a (= a;, y, z) com- 
ponent of the spatial center of the single-particle wave 
function \Lpi) is ^^7=^- When sufficiently large values 



are chosen for 



„/3 



and Wcnsti the resultant values (i and 



dc„-c„ of energy variation become /?, cic„-c„, respec- 
tively. We constrained the (iQ,_36Ar and (ii2Q_28gj values 
for the d-constraint. In each calculation of energy vari- 
ation, we constrained one of /?, d^-aeAr and di2c_2sgi. A 
detailed explanation regarding the d-constraint may be 
found in Ref. [231 . 

The energy variation with the AMD wave function 
is carried out using the frictional coohng method [1^. 
The time evolution equation for the complex parameters 
Zi , ai and jii is 



dt 



(i = l,2,... ,^), 



(8) 



where Xi is Z^, ai or and that for the real parameters 



, Vy and Vz is 



, (cr = x,y,z). 



(9) 



The quantities fix and fi^, are arbitrary positive real num- 
bers. The energy of the system decreases as time pro- 
gresses, and after a sufficient number of time steps, we 
obtain the minimum energy state. 

After the constrained energy variation for |$'^), we 
superposed the optimized wave functions employing the 
quadrupole deformation parameter p and the distances 
between the centers of mass among clusters (ic,n-c„ for 
Cm-Cn configurations as the generator coordinate. 



1$ 



M 



MK 



K 




where P^jx '^^ parity and total angular momentum 
projection (AMP) operator, and |<I>f) and |$^'^"~'^°) are 
optimized wave functions with (3- and dc„-c„"Constraints 
for the constrained values ^'^^ and _^ respectively. 
The integrals over the three Euler angles included in 
^iiK evaluated by numerical integration. The co- 



efficients f^J^ and /j^" are determined by the Hill- 
Wheeler equation. 
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M I 



M I ^Af , 



= 0. 



(11) 



C. Squared Overlap 



We defined squared overlap SO to estimate the over- 
lap of a specific model space in the superposed wave 



TABLE I: Expectation values of kinetic and potential energies 
around the ground state, ND and SD minima. K, Vnn, Vls 
and Vc denote kinetic energy, central, spin-orbit and Coulomb 
potentials. 



interaction fi 


total 


K 


Vnn 


Vls 


Vc 


Gogny DIS 0.00 


-340.9 


634.2 


-1046.2 


-0.1 


71.2 


Skyrme SLy7 0.00 


-340.7 


637.9 


-1050.0 


-0.0 


71.5 


Gogny DIS 0.39 


-326.7 


654.6 


-1033.7 


-18.4 


70.9 


Skyrme SLy7 0.39 


-326.8 


653.7 


-1031.9 


-19.4 


70.8 


Gogny DIS 0.62 


-322.2 


670.9 


-1029.9 


-33.6 


70.4 


Skyrme SLy7 0.61 


-322.3 


676.2 


-1035.9 


-33.3 


70.7 



function. Suppose the non-orthonormalized wave func- 
tions y^'f) span the functional space {AT}, for example, 

{A} = {|$f)}, {|<i>P^-)} or orthonor- 
malized wave functions |$^) are obtained by performing 
unitary transformations from I'E'j'^), 



■Jap- 



(12a) 
(12b) 



Using the l^*^), squared overlap between 1$) and the 
functional space {A}, SO is defined as 



(13) 



III. RESULTS AND DISCUSSIONS 

A. Various Structures Obtained with /3- and 
d-constraints 

We performed the energy variation after the projection 
to the positive-parity state imposing two different kinds 
of constraints, /3- and d-constraints. 

Figure [1] (a) shows the obtained energy curves as a 
function of matter quadrupole deformation (3. By apply- 
ing the /3-constraint we obtained energy curves for DIS 
(solid line) and SLyT (dotted line). Both forces give quite 
similar curves that have three local minima or shoulders 
at /3 ~ 0, 0.4 and 0.6. As shown in Table. HI they also give 
approximately the same kinetic and potential energies 
around each minimum. Therefore, we mainly discuss the 
DIS result and make some comments on the differences 
between DIS and SLy7 below. The lowest minimum at 
/?=0 corresponds to the spherical ground state and the 
two minima at /? ~ 0.4 and 0.6 correspond to the ND 
and SD states respectively. The excitation energies of 
the ND and SD minima are approximately 14 MeV and 
18 MeV. This result qualitatively agrees with the con- 
strained SHF calculation using SLy4 [l^. On the other 
hand, the energy curve in the constrained SHF-BCS cal- 
culation with SLy6 [l^l reveals different behavior. It does 
not have the ND minimum nor shoulder, but has a SD 
minimum. The excitation energy of the SD minimum 
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FIG. 1: (a): The energy curves for the positive-parity states 
obtained using j3-, d„_36Ar" s-nd di2c_2ssi-constraints. The 
energies are plotted as functions of quadrupole deformation 
/3. The solid (dotted) line shows the resuhs for DIS (SLy7). 
(b): Projection of the energy curves onto the /3-7 plane. The 
solid (dotted) line shows the results for DIS (SLy7). 



is much smaller than ours and that of Ref. [l9|. This 
difference may be due to the strong pairing correlation 
reported in Ref. H^- 

In the present calculation, we do not make any assump- 
tions regarding the spatial symmetry of the wave function 
nor put any constraints on the quadrupole deformation 7. 
7 therefore has the optimal value for each given value of 
/3. Figure [T](b) shows the projection of the energy curve 
onto the /3-7 plane. It shows that in most regions the sys- 
tem is triaxially deformed and the degree of the triaxial 
deformation greatly changes as a function of (3. Start- 
ing from the spherical ground state, the system rapidly 
changes to oblate deformation around /3=0.15-0.25 via a 
small prolate deformation. Then, it changes to triaxial 
deformation around /3=0.35 where the ND minimum ap- 
pears and the system has the largest triaxiality. With a 
further increase of /3, 7 decreases gradually from 7=30° 
at /3=0.35 to 7=15° at /3=0.70. The SD minimum that 
appears at /3=0.60 also has a large triaxial deformation 
7=15°. 





SD: zy-plane SD; 2:a;-plane SD: sy-plane 




FIG. 2: Density distributions of the ND and SD minima ob- 
tained using the /3-constraint. The crosses in the figure show 
the centroids of the wave packets. The deformation param- 
eters are (/3,7) = (0.39,25.2°) and (0.62,14.5°) for the ND 
and SD minima, respectively. 



The deformation of the (3-^ curve for SLy7 shows sim- 
ilar behavior except in the region of /3=0. 15-0.30. In the 
(3 ^ 0.30 region, the deformation obtained for SLy7 is 
triaxial, which is different from the oblate deformation 
with DIS. The wave functions in this region did not af- 
fect the ND and SD bands, nor their side bands as shown 
later. 

Figure [2] shows the density distribution of the ND and 
SD minima, and confirms their large triaxial deforma- 
tion. It also reveals the deformed mean-field nature of 
the ND minimum and implies the existence of a rela- 
tionship between the SD minimum and the cluster struc- 
ture. At the limit when all — > (all centroids of 
single particle wave packets go to the origin of the co- 
ordinate frame), the AMD wave function is identical to 
an eigenstate of the deformed harmonic oscillator whose 
oscillation number is given by the relation L0„=2hi'„ / M . 
In contrast, when the system has a cluster- like struc- 
ture, the Zi are separated into several groups to describe 
cluster subunits. In the case of the ND minimum, the 
centroids of the single particle wave packets are located 
around the origin and the wave packets are strongly de- 
formed, 



v,^ = (0.16,0.14,0.11), suggesting that 
the nature of the mean- field is deformed. In the case 
of the SD minimum, the wave packets are further de- 
formed, iv^^Vy^Vz) — (0.17,0.15,0.10), and reveal the 
aspect of a triaxially deformed mean-field. In the case 
of the SD minimum, the centroids of the single particle 
wave packets appear separated into two (28-1-12) or three 
((12-f 16)-|-12) groups and the density distribution reveals 
an octupole deformation implying that the SD minimum 
also has an asymmetric cluster-like nature. It will be 
shown using the d-constraint that this has a significant 
overlap with the ^^C-^*Si cluster structure. 

There are few studies on the triaxial deformation of 
^"Ca, but some theoretical work has suggested the triax- 
iality of the ND 0, i] and SD [ll| states. The present 
calculation with the /?-constraint has shown that most 
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of the states on the energy curve have a triaxial defor- 
mation. We therefore consider it indispensable to study 
the issue without the assumption of spatial symmetry in 
order to understand the excited states of ^"Ca. 



Q-^^Ar: type (a) a-^^Ar: type (b) '^C-^^Si 




FIG. 3: Density distributions of intrinsic states obtained with 
the d-constraint. d^_3exr is fixed to 5.0 fm for a-'^^Ar (a) and 
(b). di2c-28si is fixed to 4.0 fm for ^^C-^Si. The crosses in 
the figure show the centroids of the wave packets. 

With the /3-constraint, we do not find spatially de- 
veloped clustering at the ND and SD minima, although 
the density distribution of the SD state implies its rela- 
tionship to clustering. We have therefore applied the d- 
constraint. We discuss the results obtained with c?„-36Ar- 
and cii2Q_28gi-constraints. Other combinations of clusters 
such as 2a- ■^^Si have been also applied, but they have 
appeared at comparatively high excitation energy levels 
and were not involved with the ND and SD states even 
after the GCM calculation. 

By applying the c?Q-36Ar-constraint, we have obtained 
an excited energy curve above the energy curve obtained 
with the /3-constraint (Fig. [T|). Two different kinds of 
structure appeared on the energy curve. In both cases, 
the system has well developed a-^^Ar cluster structures 
as may be clearly seen in their density distributions (Fig. 
ini (a) and (b)). The difference between them is the 
orientation of the axis of symmetry in the oblately de- 
formed ■^^Ar cluster. The first is denoted as type (a). 
In this type, the axis of symmetry for ^^Ar is perpen- 
dicular to the vector that connects the a and ^^Ar clus- 
ters. Therefore the whole system has a triaxial defor- 
mation. This type of structure is favored under d^-seAr- 
constraints with large dc-36ArValues and appears in the 
region of (3 =0.45-0.7. The second kind of structure is 
denoted type (b). In this type, the axis of symmetry in 
the '^^Ar cluster is parallel to the vector that connects the 
a and '^^Ar clusters, resulting in an axial deformation of 
the system. This type is obtained when the inter-cluster 
distance is restricted to a shorter distance {da-sej^^—A.b- 
5.5 fm) and appears in the region of /3=0. 15-0.3. It is 
interesting that type (b) is bound deeper than type (a) 
for shorter inter-cluster distances and '^^Ar changes its 
orientation as the inter-cluster distance becomes larger. 
We assume that ^^Ar changes its orientation to make 
the overlap and potential energy between the a and ^^Ar 
clusters larger. 

By applying a (ii2c_28gi-constraint ((ii2c_28si=4. 0-6.0 
fm), we obtained a strongly deformed and excited en- 
ergy curve that appears in the region of /3=0.55-0.8. The 
system has a prominent ^^C-^*Si cluster structure and is 
triaxially deformed as shown in Fig. [3) In this case, we 



did not find a change in the orientation of the clusters. 

d-constraints have generated excited energy curves in 
which the system has prominent cluster structures. The 
fact that these wave functions are mixed with the wave 
functions obtained with the /3-constraint and play an im- 
portant role in describing highly excited bands is dis- 
cussed below. 



B. Angular Momentum Projection 
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FIG. 4: The energy curves of the J'" = 0+ and 2+ states 
projected from wave functions obtained with the /3-constraint 
with DIS. 



The wave functions obtained above are projected onto 
the eigenstate of the total angular momentum. The 
J'^—O'^ and 2+ states obtained with the /9-constraint are 
shown in Fig. [4l AMP reduces the excitation energies of 
the deformed states. For example, the states of 

ND and SD are lowered by approximately 6 and 9 MeV 
respectively. As a result, they almost degenerate in terms 
of energy. We have obtained two 2+ states for each given 
value of (3 in the deformed region by the diagonalization 
of the quantum number K . These states are denoted as 
K^=0+ and 2+ according to the dominant component of 
their wave function. The presence of the state 
is due to the triaxial deformation in the f3 > 0.3 region. 
As will be discussed below, the triaxial deformations of 
the ND and SD states leads to the presence of their side 
bands, X^=2+. 

The energy of the wave functions obtained with the 
(i-constraint are also lowered by AMP. The a-'^^Ar type 
(a), (b) and ^^C-^^Si states are lowered by approximately 
5-10, 10-15 MeV and 10 MeV respectively, a-^'^^r type 
(a) and ^^C-^*Si wave functions have the com- 
ponents because of the triaxial deformation. 
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C. GCM Calculation 

1. Energy levels and deformations 

After applying the AMP, we carried out GCM calcula- 
tions. For the GCM basis, we adopted 15 /3-constrained 
wave functions /? = 0.00-0.73, 7 c?Q_36Ar-constrained wave 
functions of type (a) with da-ssAr = 4.5-9.0 fm, 3 da-^SAr- 
constrained wave functions of type (b) with da-sBj^^ = 
4.5-5.5 fm, and 3 (ii2c_28gj-constrained wave functions 
with di2Q_28g[ — 4.0-6.0 fm. Then we obtained the fi- 
nal GCM wave functions by diagonalizing Hamiltonian 
and norm matrices for the parity and angular momen- 
tum projected states with K'^ = 0"*" and ±2+ for the 28 
independent wave functions. 

The theoretical energy levels of the GCM states and 
the experimental levels are shown in Fig. [S] The energies 
of the simple AMP for the main components in the band 
heads are also given. In most levels, the GCM states 
gain 1-2 MeV from the simple AMP. Experimentally, the 
K"" ^ 0+, 2+ and 0+ bands built on J" =0+, 2+ and 
0;^ states are known to be the ND band {K'" = ^nd)' 
side band of the ND band {K" = 2^j^) and the SD band 
{K"^ — Ogj-,). In the result of the GCM calculation, we 
obtained three K"^ = 0+ and two K"^ = 2+ bands in low- 
lying states above the ground state. We assigned the first 
and second K'^ = 0+ bands to the observed K'^ — Oj^j-, 
and Ogp bands respectively, and the first = 2+ band 
to the observed K"^ — 2^j~^ band, because the theoret- 
ical moments of inertia and electric transition strength 
B{E2) of these bands correlate well with the experimen- 
tal data for the corresponding bands as discussed below. 
The second K"^ — 2+ band in the results is regarded as 
the side band K'^ = 2^-^ of the SD band. We denote 
the third theoretical K'' = 0+ band as the K"" = 0^^ 
band because of the oblate shape. We also obtained the 
= 0+ and 2+ bands in highly excited states with 
large a-^^Ar cluster structure components, which are the 
candidates for a-^^Ar higher-nodal bands, K"^ = 0^^^ and 
2^^, observed experimentally [isj . 

Let us consider the GCM results of the low- lying states 
(ground, ND, SD and oblate states) in more detail by 
analyzing the squared overlap between J'^ = 0+ states 
and the AMP states from the /3-constrained wave func- 
tions. The ND state's K"^ — Ojjjp and 2^j^ bands are 
mainly constructed from the /3-constrained wave func- 
tions around the ND minimum occurring in the (3 ^ 0.4 
region. Their band-head states have a maximum overlap 
of approximately 85% for the triaxially deformed state 
with {/3, 7) = (0.39,25.2°). There is no mixing between 
the ND state and the ground state. The squared over- 
lap with the /3-constrained wave functions is almost un- 
changed up to high spin states along the K'^ = Oj^j-, band. 
On the other hand, in the K"^ = 2^^) band, the squared 
overlap of each component changes with the increase in 
spin, which implies the change of structure in this band. 
In the result of the GCM calculation, the theoretical ex- 



TABLE II: The squared overlap {SO) of /3-constrained wave 
functions and d-constrained wave functions. 
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0.61 
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a-^'^Ar 




0+ 


0.41 


0.49 a-^^'Ar 


higher-nodal 
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0.45 


0.50 






4+ 


0.52 


0.50 






2+ 


0.38 


0.55 






3+ 


0.52 


0.45 






4+ 


0.59 


0.39 



citation energies of the band heads in the K"^ = OjJ^j-, 
and 2'^Y) bands are 11.2 MeV and 12.1 MeV respectively. 
These are much higher than the experimental excitation 
energies. However, the rotational energies in the band 
are reproduced well, as shown in Fig. [6l 

The SD states, K"^ = Ogj-, and 2gj3, are constructed 
mainly by /3-constrained wave functions around the 
SD local minimum. The main component is the /3- 
constrained wave function with triaxial shape, (/3, 7) — 
(0.62,14.5°). The squared overlap is more than 90% in 
the band-head states, and almost unchanged up to high 
spin states along the K"^ = Og^ and 2gj-| bands. A small 
degree of mixing between the ND and SD bands is seen 
in the f3 ~ 0.5 region. The result of the GCM calcu- 
lation reveals the band-head energies of the = O^D 
and 2gQ bands are 12.1 MeV and 14.0 MeV respectively. 
Although the present calculations overestimate the ex- 
perimental excitation energies of the K'' = 0+^ band as 
well as the ND bands, the rotational energies in the band 
are reproduced well as shown in Fig. [S] 

The oblate bands are dominated by the /3-constrained 
wave functions in the /3 ^ 0.3 region. The main 
component is the oblate wave function with (/3, 7) = 
(0.33,51.2°). The oblate band {K'' = 0+^) has not been 
discovered experimentally yet. Our prediction of the 
oblate band is consistent with spherical-basis AMD [2l| 
as well as SHF calculations ^] . 



2. Cluster components 

In this section we discuss the contribution of the cluster 
wave functions in the ND and SD states. In Table HIl we 
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FIG. 5: Energy levels in ^°Ca. The left hand side represents experimental values, and the right hand side theoretical values. 
The energies of superposed wave functions and that of the main component are represented. 




angular momentum / 



FIG. 6: The rotational energies E^ot = E — Eq+ of the ND 
and SD bands as functions of the angular momentum J( J+1). 



list the squared overlap {SO) values of the GCM states 
with the a-^^Ar configuration space given by the set of 
da-seAr-constrained wave functions, and those with the 
i2(-;_28gj configuration space, as well as the SO values for 
the model space of the /3-constrained wave functions. The 
definition of the specific functional space is explained in 
mi CI The squared overlap between the ND state and f3- 



constrained wave functions is 99%. This means that the 
ND state is practically represented by the /3-constrained 
wave functions alone. On the other hand, it is surpris- 
ing that the ND states have a significant overlap of ap- 
proximately 40% with the da-seAr-constrained wave func- 
tions as well. This indicates that the /^-constrained wave 
functions for the ND states include the a-^^Ar compo- 
nent, though the spatially developed cluster structure is 
not seen in the density distributions. The a-^^Ar cluster 
component in the ND band is mainly consistent with type 
(a) wave functions, while the squared overlap of type (b) 
wave functions in ND states is almost negligible. The re- 
sult, that the ND state contains an a-'^^Ar cluster struc- 
ture component, is associated with the results of the a- 
36 Ar potential model calculation ^^,^^Ar{'^U,dy°Ca 
reaction [ll[ll| and a-^^Ar OCM calculation 

The SD states are dominated by the /3-constrained 
wave function as well as the ND states. However, the 
SD states also have a large overlap with the c?i2(;;_2ssi- 
constrained wave functions as reflected by the SO of ap- 
proximately 60%. We found that the /3-constrained AMD 
wave functions for the SD states also include the cluster 
components, even though the cluster structure is not vis- 
ible in the density distributions. The ""^^C-^^Si cluster 
configurations also make a significant contribution to the 
energy of the SD states. In particular, the energies of 
the SD states gain 1-2 MeV due to the mixing of ^^C- 
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^®Si cluster structure wave functions. This is associated 
with the resuhs of spherical-basis AMD [2l| according 
to which the SD states have a ^^C-^*Si cluster structure 
configuration. 



3. a-^^ Ar higher-nodal states, K'" = Q'^^ and 2^^"^^ hands 

As shown in Fig. [5l we obtained K"^ = 0+ and 2+ 
bands with large a-'^^Ar cluster components in the exci- 
tation energy region approximately 10 MeV higher than 
the ND band. The main component of the higher a-^^Ax 
bands, K'" = Oj|„ and 2+^, are c?Q_36Ar-constrained wave 
functions obtained with a large distance da-^exr = 6.0 fm. 
We assume that this corresponds to the a-'^^Ar higher- 
nodal band observed with the '^^Ar(^Li,d)^°Ca reaction 
[l5l |. where the fragments of the J'" — 0+ state for this 
band were reported around 8 MeV above the band-head 
of the ND band. 



0.40 




4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 

4-»Ar [fm] 



FIG. 7: The squared overlap of the J" — 0^ states in the 
ND band [K"^ = OjJ^p) and the a-^^Ar higher-nodal band 
{K'^ — Ojljj) with d^_36Ar-constrained wave functions. 



The squared overlap between dQ,-36Ar-constrained wave 
functions and the J'^ = 0+ states in K'^ — O^p and OJ^j^ 
bands are shown in Fig. [71 The squared overlap for is 
suppressed in the small d^-sej/^j. region and has a peak at 
da-sBAr ~ 6 fm, which demonstrates the nodal property of 
the inter-cluster motion in this band. It can be assumed 
that the a-^^Ar higher-nodal states arise from the inter- 
cluster excitation built on the ND states. In other words, 
the significant component of the a-^^Ar cluster structure 
in the ND states must be essential for the formation of the 
higher-nodal states, because the cluster component in the 
small distance region is possible if the ND states contain 
no cluster component. This means that the higher- nodal 
states appear as a consequence of the orthogonality to 
the cluster components with a small distance d, which 
are already contained in the ND states. This situation is 
similar to the relationship between the ground band and 



TABLE III: Theoretical and experimental B{E2) of ''"Ca. 
The units of transitions are Weisskopf units, B{E2)w.u. = 
8.12 e^fm*. The asterisk-marked experimental B{E2) val- 
ues are taken from Ref. [Tsl . and other values are taken from 
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2.6 ±0.7 






4+ 2+ 


< 0.05 


2.6 ±0.6 




2+ 
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0.2 


22 ±6 



Q!-'*°Ca higher-nodal band in Ti suggested by Kimura 
et al. with deformed-basis AMD [28l |. 



4- Eleetrie transitions 

Here, we investigate the electric quadrupole strengths 
and discuss the band structure. The theoretical and ex- 
perimental values of E2 transition strengths B{E2) are 
shown in Fig. [51 The B{E2) values for intra-band transi- 
tions are remarkably strong, and the transitions between 
the K"^ = O^Y) and = 2^^) bands, and those be- 
tween the K'^ = OgQ and 2gp bands are also strong. 
This reflects the side-band features of the = 2^^ 
and K'^ = 2^-^ bands, corresponding to K"^ = O^j-, and 
K'^ ~ 0gj3 respectively. The inter-band transitions be- 
tween ND and SD states and those from ND or SD states 
to the ground state are underestimated. This is because 
mixing of wave functions among these bands is small in 
our calculation. On the other hand, we predict the tran- 
sitions from the oblate — O^j, band to the ND bands, 
K'' = 0+J3 and 

2nd' because of the mixing among these 

bands. 

Detailed comparisons between the theoretical values 
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FIG. 8: Experimental and theoretical B{E2) values. Units are Weisskopf units, B{E2)w.u. = 8.12 e^fm'* for ''^Ca. In 
theoretical values, transitions stronger than 4 4 B{E2)w. u. are presented. The experimental data are taken from Refs.[l3 and 
[2^ . "unknown" indicates the transition has been observed but B{E2) has not been measured. 



and experimental data of the B{E2) values are shown 
in Table [nil The B{E2) values for the intra-band tran- 
sitions in the ND and SD states are reproduced well, 
except for the 6+ 4+ transition when = O^j-,. The 
small value of the experimental B{E2; 6+ — > 4+) is over- 
estimated by the calculation, although the error in the 
experimental value is rather large. 

Quadrupole moment 

The theoretical and experimental values of quadrupole 
moments deduced from the E2 strengths for the K"^ = 
^ND ' ^^"^ ^SD bands are listed in Table IIVI The 

results of other theoretical studies are also shown in the 
table. The intrinsic quadrupole moment Q is defined by 
the B{E2) value as follows. 



Q 



'16tt B{E2, I, If) 
5e2" {I,K20\IfK)^ ' 



(14) 



where {IiK20\IfK) is the Clebsh-Gordan coefficient. 

The theoretical quadrupole moments of the ND and SD 
states are around 130 fm^ and 220 fm^ respectively. The 
values are almost unchanged up to high spin. This means 
that the ND and SD bands consist of approximately rigid 
rotor states. These results are consistent with the rigid 
rotor-like spin dependence of the excitation energies of 
the ND and SD states shown in Fig. [Sj The theoretical 
Q moments are consistent with experimental values eval- 
uated from life-time measurements and branching ratios 
(I) except for the 6+ 4+ transition in the ND state. 
The small value of the experimental Q moment for the 



6+ — > 4+ transition might suggest some structure change 
in high spin states in the ND band. 

We compared our results with those of other theoret- 
ical studies. The Q values in (4) a-^^Ar OCM are con- 
sistent with the present values for the ND states. This 
may indicate that the deformation of the ND states in 
our calculations is similar to that calculated with the a- 
•^^Ar cluster model. In (1) the spherical-basis AMD and 
(3) the shell model, smaller Q values for the K"^ = 0+ 
states were obtained by comparison to our model. In (2) 
SHF-BCS -I- GCM, the Q values are also smaller than our 
results for = 0+ states. In particular, the Q moments 
for the 4+ 2+ transition in both the ND and SD states 
are remarkably small. This shows that the ND and SD 
bands obtained in the SHF-BCS -I- GCM calculation are 
not rigid rotor-like. This seems to be inconsistent with 
experimental results regarding the rigid rotor-like prop- 
erty, which has been found in rotational energies and E2 
transitions. In recent experimental observations, the Q 
values have been obtained by globally fitting all transi- 
tions. They ((Il)-a and b) are smaller than the values of 
the life-time measurement (I) and also smaller than the 
present results. Further analysis requires more detailed 
measurements of the E2 transition strengths. 

Triaxial rotor 

In this section we discuss the characteristics of the side 
bands, K'' = 2+, in the ND and SD states. In a naive 
collective model, a K"^ = 2+ band can be described by 
either a triaxial rotor or a 7 vibration. In order to identify 
the mode of the side bands, it is useful to analyze the 
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TABLE IV: Intrinsic quadrupole moments of the ND band and the SD band extracted from B(E2). For theoretical study, (1), 
(2), (3) and (4) are taken from the results of spherical-basis AMD HFB-BCS + GCM ;20], shell model ;22] and a-^^'Ar 
OCM [13j, respectively. For experimental data, (I) is obtained from life-time and branching ratio, asterisk-marked values are 
taken from Ref. [l^ and others are from Ref. [2^. (Il)-a and (Il)-b are taken from Refs. [3 and|23 respectively, which are obtained 
through a global fitting of all the transitions. 
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ratio of E2 strength for K'^ — 
B{E2]{2+,K'^ = 2+) 



R = 



2+ to K"" = 0+ transitions 
^ (0+,^'^ =0+)) 



B{E2; {2+,K'' ^ 2+) ^ (2+, K'' = 0+)) ' 



(15) 



where B{E2;{Jf .R-^) {JJ,KJ)) is the value of 
quadrupole electric transition strength from a state 
in the Kf band to a JJ state in the KJ band. The R 
value is set to the constant value 7/10 in the 7 vibration 
mode, while it is a function of the static 7 parameter in 
the triaxial rotor mode [29] , 



ff°'(7) 
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20 



3-2sin2(37) 



9-8sin2(37) 



sin2(37) 
9- 8sin2(37)' 



(16) 

The function i?''°*(7) is plotted as a function of 7 value 
in Fig. [9l The value of R'^°^ is 0.7 for the axial symmetric 
limit, 7^0° and 60°. As the triaxiality increases, R""^ 
decreases, going to zero at the peak of triaxiality, 7 = 
30°. Regarding the ND state, the theoretical value R 
of the E2 strength ratio is nearly equal to 0, and the 
experimental value R — 0.06 is also close to 0. Both the 
values agree well with i?''°'(7 ~ 30°) but are far from that 
of the 7 vibration limit, 0.7. In the present results, the 
7 value of the main component of the ND state is 25.2°, 
which is consistent with 7 30° for the triaxial rotor 
derived from the experimental R value. We therefore 
conclude that the mode of the observed K"^ = 2^^ band 
is a triaxial rotor, not a 7 vibration. 

Regarding the SD state, the theoretical value of R 
is equal to 0.51, and the 7 value of the main compo- 
nent of the theoretical SD state is 14.5°. These values, 
(7,i?) = (14.5°, 0.51) can be described by the relation 
_R''°'(7) in Eq. for the triaxial rotor. This indicates 
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FIG. 9: The values of the ratio of E2 strength _R as a func- 
tion of 7 value. The dot-dashed line and the solid line show 
the limiting values of the 7-vibration and the triaxial rotor 
respectively. The i?ND(~ 0) and -Rsd(= 0.51) values are the 
theoretical R values for the ND and SD states derived from 
the E2 transition strengths respectively, while that of 

the experiment for the ND state. The deformation parame- 
ters 7nd(= 25.2°) and 7sd(= 14.5°) for the main components 
of the ND and SD states are also shown. 



that the case of K'^ = 2gp in the SD state is the side 

The K'' = 2^0 band has not 



band of the triaxial rotor 
yet been observed experimentally. In order to confirm the 
triaxiality of the SD band, observations of the K"^ = 2^-^ 
band and the E2 strength for the inter-band transitions 
are required. 
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IV. SUMMARY 

We investigated the ground state and excited states 
of ^°Ca in the framework of deformed-basis AMD fo- 
cusing on the aspects of triaxiaUty and clustering in 
the deformed states. Superposing mean-field-type and 
cluster structure-type wave functions obtained with /?- 
and (i-constraints respectively, wc obtained the ground 
states, ND, SD, oblate and a-^^Ar higher-nodal states. 
We found both the ND and SD bands are constructed 
from triaxially deformed shapes, and we obtained the 
= 2+ side bands of these bands. The theoreti- 
cal B{E2) and moments of inertia are consistent with 
experimental data, although the excitation energies are 
higher than empirical values. By analyzing quadrupole 
moments calculated from the B{E2) and the E2 strength 
ratio R values, we found that the ND and SD states can 
be understood by a relatively rigid rotor conception of 
triaxial deformations. The ND band and it's side band 
contain approximately 40% of the a-^^Ai cluster struc- 
ture component, and the SD band and it's side band con- 



tain approximately 60% of the ^^C-^^Si cluster structure 
component. The a-^^Ar higher-nodal band was obtained 
due to the excitation of inter-cluster motion between the 
a and ^^Ar clusters. The present results suggested that 
cluster correlation will be important for deformation and 
excitation even in medium- and heavy-weight nuclei. 
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